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We study decoherence in a simple quantum mechanical model using two approaches. Firstly, 
we follow the conventional approach to decoherence where one is interested in solving the reduced 
density matrix from the perturbative master equation. Secondly, we consider our novel correlator 
approach to decoherence where entropy is generated by neglecting observationally inaccessible cor- 
(-H , relators. We show that both methods can accurately predict decoherence time scales. However, 

^ • the perturbative master equation generically suffers from instabilities which prevents us to reliably 

calculate the system's total entropy increase. We also discuss the relevance of the results in our 
, quantum mechanical model for interacting field theories. 

^ ■ PACS numbers: 03.65.Yz, 03.67.-a, 98.80.-k,03.65.-w,03.67.Mn 

I. INTRODUCTION 
A. The Correlator Approach to Decoherence 

o ■ 

\ Recently, we proposed a novel approach to decoherence based on neglecting observationally inaccessible correlators 
[Hi 13 J ^1^0 see 0. Just as in the conventional approach to decoherence [j-Q, we split our world into a distinct 
system, environment, and observer. Interactions between system and environment in quantum field theory generate 
LiJ I non-Gaussian correlators. However, from the observer's perspective, these correlators are generically small in the 
' perturbative sense. The observer's inability to measure the information stored in these higher order correlators 
generates a certain amount of entropy which turns our quantum system into a classical stochastic system. We applied 
our "correlator approach" to decoherence to an interacting quantum field theoretical model with an interaction 
term of the form h(j){x)x^ {x) , by solving for the renormalised statistical propagator in an out-of-cquilibrium setting. 
; ^ Giraud and Serreau [3| performed a similar analysis in a A0^(a;) field theory. Starting point is the well known von 
I Neumann entropy: 

S'vN = -Tr[/51n(/5)] , (1) 
where p denotes the density operator which in the Schrodingcr picture satisfies the von Neumann equation; 

th^^p^[H,p]. (2) 

By restricting ourselves to a Gaussian density matrix (and thus neglecting all non-Gaussian information that is 
contained in the full density matrix), one can easily relate the Gaussian von Neumann entropy to the statistical 
propagator of our field theoretical system [2],Q . This procedure can be improved by including non-Gaussian corrections 
to the entropy Q • Indeed, the Gaussian von Neumann entropy is in most cases no longer conserved, unlike the total 
von Neumann entropy in equation ([T]). 

We thus argue to solve directly for the Gaussian correlators characterising the Gaussian properties of our system. 
A sophisticated field theoretical framework exists to renormalise the statistical propa gato r, include perturbative 
corrections and study genuine time evolution in an out-of-cquilibrium setting (see e.g. (ol Il0j|). 
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B. The Master Equation Approach to Decoherence and its Shortcomings 

111 the conventional approach to decoherence one is interested in solving for the reduced density matrix which is 
obtained by tracing over the environmental degrees of freedom in the full density matrix: 

PrcdW = Tr£[p(t)]. (3) 

The unitary von Neumann equation ^ transforms to a non- unitary master equation for prod- The master equation 
is in principle equivalent to the von Neumann equation if no approximations arc made. This is for example clear 
from the influence functional method. In other words: the exact master equation is just as hard to solve as the von 
Neumann equation in interacting quantum field theories. Therefore, one usually relies on perturbative methods to 
simplify the exact master equation and to derive a perturbative master equation. 

The conventional approach to decoherence suffers from both theoretical and practical shortcomings. From a the- 
oretical point of view, it is disturbing that pi-cd is obtained from a non-unitary equation despite of the fact that the 
underlying theory is unitary^. Practically speaking, the (exact or perturbative) master equation is still so hard to 
solve that even the most basic field theoretical questions have never been properly addressed: no well established 
framework exists to include perturbative corrections to a reduced density matrix (see however the pioneering work in 
quantum mechanical cases pJHlSl l) nor has any reduced density matrix ever been renormalised^ . 



C. A Clean Comparison Between the Two Approaches 



As a solution to the perturbative master equation in a field theoretical setting has so far not been derived, let 
us study a simple quantum mechanical model. Quantum mechanics thus provides us with the ideal playing field to 
compare the two approaches to decoherence. In this paper we study a well known textbook example |13l . Il5l4l7l | of 

-I- 1 simple harmonic oscillators. Here, one oscillator, x, plays the role of the system and the other N oscillators, 
Qn, 1 < n < N , play the role of the environment. The oscillators will be coupled quadratically: A„a;g„. In quantum 
mechanics we do not have to worry about divergences and moreover, since we have chosen our model to be Gaussian, 
we can actually compute the reduced density matrix by tracing over the environmental degrees of freedom. Indeed, all 
the practical shortcomings to the conventional approach to decoherence as listed above have thus been circumvented. 

Furthermore, it is important to realise that the entropy in our correlator approach to decoherence is not generated by 
neglecting non-Gaussianities as our quantum mechanical model is Gaussian. Instead, the correlation entropy between 
the system and environment that builds up due to the coupling between the two is neglected. The correlation entropy 
in this paper thus plays the role of the non-Gaussian contributions to the entropy in a proper field theory. In this 
paper we present two calculations: 

• We calculate the evolution of the entropy by solving for the reduced density matrix from the perturbative master 
equation; 

• We calculate the evolution of the entropy by solving for the Gaussian correlators directly from which we calculate 
the Gaussian von Neumann entropy. 

In particular we identify regions in parameter space where the perturbative master equation performs well and where 
it does not. In section [IT] we introduce the model and in section Hill we outline two derivations of the perturbative 
master equation. In sections HVl and IVl we analyse the evolution of the entropy for = 1 and A^ = 50, respectively. 

II. THE MODEL 

A. BiHnearly Coupled Simple Harmonic Oscillators 

We consider the following total quantum mechanical action: 

S [x, {g„}] = / dtL [i, {g„}] = / d< {Lg [i] + Le + Li„t [i, , (4) 



^ The energy is not conserved as a consequence. This means that this basic method to check a particular numerical evolution for Prcd is 
not available. 

^ For example in [l4j the decoherence of inflationary primordial fluctuations is investigated using the master equation approach however 
renormalisation is not addressed. 
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where: 



1 ~2 2-2 

-mx mujQX 



Le [{qn}] = ( " I 



n=l 



N 



Lint [x, {q„}] 



Here i = denotes the system coordinate and g„ = qn{t) with n = 1, ■ 
is useful to perform the rescaling: 



after which the Lagrangian in equation Q becomes: 



mx 



(5a) 
(5b) 

(5c) 

, N are the environment coordinates. It 

(6a) 
(6b) 



L[x,{qn}] = -i^ 



^22 



^ /I 1 \ ^ 



71=1 



(7) 



n=l 



where A„ = A„/^77im„. Throughout the paper we will set (x) = and (qn) ^ 0, 1 < n < N . In order to 
understand how to impose initial conditions later on, it is useful to explicitly derive the statistical propagator in 
thermal equilibrium. In thermal equilibrium the Schwinger-Keldysh propagators for e.g. the coordinate x are given 
by: 



iA++{t;t') 
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(8b) 
(8c) 
(8d) 



where /3 = l/{kBT), with fc^ the Stcfan-Boltzmann constant and where we defined At = t — t' . The causal and 
statistical propagators now follow as the sum and the difference of the two Wightman propagators: 



A%t;t') = A-+(i;t')- A+-(i;t') = - 
F{t;t') = ^{iA-+{t;t') + tA+-{t;t') 



sin(cjoAt) 

cos(a;oAt) 
2uJo 



coth 



/3ujo 



(9a) 
(9b) 



B. The Von Neumann Entropy in Quantum Mechanics 

1. A = 

It is important to realise that since the Lagrangian in equation ([7]) is Gaussian, no non-Gaussianities are generated 
by the coupling A. For a free quantum mechanical system (A = 0), it is simple to find the total von Neumann entropy 
in equation ([T]): 

•SvN = -Tr[pf logpf] , 

where the subscript "f" for "free" reminds us that all interactions or couplings are switched off. The density operator 
of a quantum mechanical Gaussian state centered at the origin is of the form: 



Piit) 



oo J — oo 



d.T / Ay\x)pi{x,y\t){vi\ 



(10a) 
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where: 



Pf{x, y; t) = M{t) cxp \-a{t)x^ - h{t)y^ + 2c(i)xy] , 



(10b) 



and where a = a(i), h = b{t) and c = c{t) are determmed from the von Neumann equation. Moreover, from = pf 
it follows that b* ~ a and c* = c. When c 7^ the density matrix is mixed and entangled. The normalisation J\f is 
obtained from requiring Tr[/5t] = 1: 



2(aR - c) 



(11) 



provided that c < or where = 5ft[a]. We can furthermore derive the three non-trivial Gaussian correlators 
characterising our system: 



/oo 
-00 



x\pfx'^\x) 



1 



4(aR - c) 



ai 



2(aR - c) 



Or - C 

We used (xIpItA) = —idxixltp). These relations can be inverted: 



ai = 



an 



2(.t2 
A2 + 1 



We can straightforwardly find the von Neumann entropy by making use of the replica trick (see e.g. 0): 



DvN = — 7: — in 



A- 1 



In 



A- 1 



where: 



A^ - 



Qr + c 
Or - c 



[x'){f 



1 \ 2 



4 [F{t;t')dtdt,F{t;t') - {dtFit;t')f 



(12a) 
(12b) 
(12c) 

(13a) 
(13b) 
(13c) 

(14) 
(15) 



The physical meaning of A/2 is the phase space area occupied by a Gaussian state centered at the origin in units of 
h. For a pure state we have A = 1, whereas for a mixed state A > 1. The phase space area is an extremely important 
function, as it coincides precisely with the area the state occupies in Wigner space (H and is moreover conserved by 
the evolution d/di[A2] = in free theories. It is also important to mention that we can associate a statistical particle 
number to A given byn=(A — l)/2 such that the von Neumann entropy reduces to the well-known entropy of n 
free Bose particles per quantum state. 



2. A / 

As in all unitary systems, the von Neumann entropy following from equation ([1]) is conserved: 

S'vN = const . (16) 

The von Neumann entropy for the model ([7|) contains both information about the system, the environment and the 
correlations between system and environment that are generated due to the coupling: 

^vN = 5total = Ss{t) + SE{t) + SsEit) • (17) 
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where Ss and Se represent the Gaussian von Neumann entropy contained in the system and environment separately, 
and where Sse denotes the Gaussian correlation entropy. Adding the word "Gaussian" to these entropies might at 
the moment seem redundant, because our total system and environment arc Gaussian so no non-Gaussianities can 
be generated as we mentioned before. We nevertheless insist on this nomenclature in order to easily compare it with 
interacting field theories later. For a free system S that is not coupled to or in interaction with an environment E 
we have derived the von Neumann entropy in equation (|14|) and (jisp. Likewise, we can associate a Gaussian von 
Neumann entropy to our system S in the presence of a coupling A too: 

where the phase space area is given in equation psp and follows from measuring the three non-trivial equal time 
Gaussian correlators characterising the properties of the system: (i^), (p^) and l/2{{x,Px})- Note that for the 
coupled system Ss(t) is generally not conserved. A similar argument applies to the environmental oscillators and we 
can likewise associate an environmental entropy to all qn if we measure its properties. Since the total von Neumann 
entropy is conserved in equation (jl7|) . we can then find the time evolution of the correlation entropy by solving: 

SsE{t) = 5vN - Ss{t) - Ssit) = 55(0) + Se{0) - Ss{t) - SE{t) , (19) 

where we assumed that Ss b (0) = initially. What we would refer to as "Gaussian von Neumann entropy" is referred 
to as "correlation entropy" in [l8| . where they prove an iJ-thcorcm for a quantum mechanical 0{N) model. 

The reduced density matrix captures the average effect from the environmental oscillators on our system oscillator 
and is defined in equation ([3]): 

Prcd(i) = TTE[p{t)] . 

In order to quantify the amount of decoherence that has taken place, one is interested in the "reduced von Neumann 
entropy" which we can define by'^: 

Sl'^^^it) - -Tr[p,ed(i) log A-cd(i)] . (20) 



The Gaussian von Neumann entropy we defined in equation (jlSp is identical to the "reduced von Neumann entropy" 
defined above: 

55(<) = ^rN(i)- (21) 

The argument is trivial: 

/oc 
dxdq{xq\px'^\qx) = Tvslprcdx'^] , (22) 
-00 

and similar expressions hold for {{x,px}) and {p^). Since the value of all three non-trivial Gaussian correlators docs 
not depend on whether we evaluate the expectation value with the full density matrix or with the reduced density 
matrix, the value of the phase space area A(<) or indeed the entropy is identical in both cases. We show this explicitly 
using our coupled system of simple harmonic oscillators in appendix [BJ Of course, when proper (non-Gaussian) 
interactions are considered, equation pi|) is only true for the "Gaussian reduced von Neumann entropy" S^j^'°'^{t) 
based on including Gaussian correlators only. 

III. THE MASTER EQUATION 

In this section we outline for completeness two derivations of the master equation. Paz and Zurek nicely derive 
the master equation by means of pcrturbative methods [l6l | before they review the "infiucncc functional method" . 
The latter was pioneered by Caldeira and Leggett [l^, see also the early review of Grabert, Schramm and Ingold 



^ This entropy is not to bo confused with the entanglement entropy of for example black holes which is generated by the observer's 
inability to access the information in all of the space-time. 
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[isj , which in principle ahows us to derive a non-perturbative master equation [lT|, [20l - [2^ . In the perturbative 
approximation both derivations yield the same master equation. 

One defines the reduced density matrix by tracing over the environmental degrees of freedom in the full density 
matrix as in equation In most of the cases of physical interest, it is impossible to solve for the von Neumann 
equation of the density matrix exactly. Hence, it is impossible to perform the trace to find the reduced density 
matrix exactly. Consequently, in order to write down an equation of motion for prcd that is usually referred to as a 
"master equation" , one makes certain approximations to the von Neumann equation. For example, one can make a 
perturbative approximation or neglect the backreaction from the system on the environment. The master equation is 
non-unitary in nature which is precisely what generates the entropy. 

Let us state again explicitly at this point that there is nothing wrong with tracing over environmental degrees of 
freedom per se, it is the perturbative master equation that is used to obtain the "reduced density matrix" that fails 
to capture the correct physics at late times as we will see later. 

Let us begin by following the perturbative approach as outlined by Paz and Zurek in e.g. [l^. The Hamiltonian 
following from ([7]) is given by: 

11 ^ /I 1 \ 

H = Hs + He + Hint = ^Pl + 2^0^' + ( + 2'^"'^" ^ ) ■ ^^^^ 

n=l ^ ^ 

The perturbative solution to the von Neumann equation ([2]) is the Dyson series which is truncated at second order. 
We can now trace over the environmental degrees of freedom to obtain: 

/5LdW = ^Trs[i?L(i), /5(0)] ^ ^ dtiTr^[ffL(0, [Ant(ii), P(0)]] . (24) 

Here, the superscript I denotes the interaction picture such that Hj^^t{t) ~ UQHint{0)Ua, with Uo ~ cxp[~i/h{Hs + 
He)]- If one now assumes that initially the total density matrix is not entangled, i.e.: p{0) = psifi) ^ /5_e(0), one can 
rewrite equation (j24p in order to express it solely in terms of pl^^{t): 

Ldit) - ^Ti-E[Hltit),pldit) ® PEiO)] - ^ ^ dhTrEiMntit), [HL{ti),plodit) ® PEiO)]] (25) 

+ ^J^ dhTvE {[HUt),^i-E {[Hltiti),pldit)(^ PEiO)])® PEiO)]) . 

For the model under consideration ([23|), the master equation above reduces in the Schrodinger picture to (fi, = 1 
again) : 

1 ~ /■* 

Acd(i) = -[ffs(i),/5rcd(<)] - / dhiyiti)[x,[xi-h),p,,dit)]]-iviti)[x,{xi-h),p,,dit)}]. (26) 
« Jo 

Here, note that i(— t) = xcosiujot) — px/ujQs'miujQt) due to changing back to the Schrodinger picture from the 
interaction picture. The noise and dissipation kernels i^it) and 77 (i) read at the lowest order in perturbation theory: 

Ht) - (a(t),^„(0)})^f:A^„H^^eoth(^) (27a) 

n=l n=l ^ ^ 

nit) - E^([^»(^)'^"(0)])-Et^^^- (27b) 

71—1 n—1 

One thus finds: 

/5rcd(i) = -i[Hsit) + ^n^it)x^,prcdit)] - i-fit)[x,{p,,p,,dit)}] - Dit)[x, [a;,prcd(t)]] - fit)[x, [p„/5,ed(i)]], (28) 
where the frequency "renormalisation" flit), the damping coefficient 7(t) and the two diffusion coefficients Dit) and 
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f(t) are given by: 

N 

y 



/■* ^ -A^ 

—2 / dii77(ii) cos(aJoii) = / — 2 " — 2T t"^" (cos(a;oi) cos(w„t) — 1) + sin(a;ot) sin(a;„i)] (29a) 

Jo ^n[i^o ~ ^n) 

7(t) = [ dtirjiti)^^^^^—^^ ^ 'y ^ 7TT [i^n cos{uint) sin{u!ot) - ujo cos{ujot) siii{u}nt)] (29b) 

Jo ^^0 2wot^„(a;^ -w^) 

D(t) ^ / diii^(^i) cos(cjoii) = > . „ — [wq cos(cj„i) sin(a;ot) — w„ cos(cjoi) sin(w„i)] (29c) 

Jo f-; 2a;„(a;^-w2) 



A2 coth 



f{t) ^ ~ [ dtiiy{ti f^^^^°^'^^ =^ V r 2 ^ N (cos(woi) cos(cj„t) - 1) + w„ sin(woi) sin(cj„t)] . (29d) 

Jo ^0 2wot^«(w^ - <) 

The crucial point is that the non-unitarity present in the master equation in the terms "f{t), D[t) and /(t) in equation 
(|28p brings about the entropy increase. This entropy increase, in turn, is argued to be a good quantitative measure of 
the process of decoherence. Surprisingly, we are not aware of any direct calculation of the entropy in this approach in 
generality. Furthermore, it is important to note that the effective coupling constant is not just given by A but rather 
by: 

This can be readily appreciated from examining the eigenfrequencies for = 1 in equation (jA4p for small couplings. 
We will return to this observation in subsequent sections, but let us at this stage already point out that when the 
system and environmental frequencies are close together, one risks leaving the regime where the master equation is 
valid. This is disturbing. When both frequencies are roughly equal one is in the so-called resonant regime, i.e.: the 
regime where the interaction between both oscillators is most effective. There is in principle nothing non-perturbative 
about the resonant regime. For example, one can study the case A'^ = 1 where wq — wi and with both A/wq ^ 1 and 
\/uj\ ^ 1. Despite of this, the coefficients in the master equation tend to become large in this regime. 

Before we solve the master equation (|28|) . let us outline the rather well known "influence functional" derivation of 
the master equation which gives the same perturbative result [H, [l^ . The full density matrix is first projected on 
position basis: 

p(x,g;2/,r,t) = / dx'dyWdr'(a;,g|e-'^*|x',q')(a;',g'l/5(0)|2/',r')(y',r'|e^^%,r) . (31) 



One can then define a quantum mechanical kernel |23| to be the position space projection of the evolution operator, 
which can be evaluated using a path integral approach: 

K{x, q, t; x\ q', 0) = (.t, q\e-'"*\x' , q') = j VxVq e'^["'«l (32a) 
K*{y, r, t- 2/', r', 0) = (y', r'|e'^*|y, r) = / VyVr e-'^^^'^^ , (32b) 



where the action follows from equation (j36p . The boundary conditions for the path integrals are given by: x{t) = x, 
x(0) = x', q{t) = q, q{0) = q' and likewise y{t) = y, y{0) = y', r{t) = r, r(0) = r' . Now, we can trace over 
the environmental degrees of freedom as usual to find the reduced density matrix. Assuming that the system and 
environment are not entangled initially, i.e.: /)(0) = /5s(0) ® pe{0), one can thus rearrange the integrals to derive: 



p,cd{x,y,t) = J dq{x,q\p{t)\y,q) ^ J dx'dy' J{x,y,t; x' ,y' ,0)ps{x' ,y' ,0) , (33) 

where: 

Jix, y, t- x\ y\ 0) = / VxVy e^^-[-l-'^«[yl^(x, y) , (34) 



and where the influence functional T is given by: 

J-(a;, 2/) = y" dQd<?'drV£;(q', r', 0) ^ VqVr e»Si=M-«SB[r]+,Si„t[x,9]-,Si„t[a,r-] (35) 

Caldeira and Leggett then derive an equation for /9ied(a;, y, i), i.e. the master equation, by carefully studying infinites- 
imal time translations /Orod(a^, 2/i ^ + 6)5 which they expand in e. 
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IV. TWO COUPLED OSCILLATORS 

We can specialise to the case of two bilinearly coupled simple harmonic oscillators by setting iV = 1 in equation ([7]) . 
The total Lagrangian we consider in this section is given by: 

This is the simplest quantum mechanical system one can think of that captures some of the physics that is relevant 
for decoherence. In this section, wc perform two computations: 

• Decoherence in the conventional approach: wc can straightforwardly solve for the master equation in this case; 

• Decoherence in our correlator approach: we can solve for the statistical propagator of the system oscillator and 
the environment oscillator by exact analytic methods from which we can extract the exact Gaussian entropy. 



The bottom line is clear: the simple case p6p allows for a clean comparison of the two approaches to decoherence. 

In appendix |B]wc compute the reduced density matrix by tracing over the environmental degrees of freedom exactly 
as the Lagrangian (|36[) is Gaussian. As we can also solve the von Neumann equation for the full density matrix, we 
can study the (unitary) time evolution of the reduced density matrix and confirm that the equality derived in equation 
(PT|) is valid. 

A. The Master Equation Approach to Decoherence 



We are now ready to evaluate the master equation (|28p. Let us project this operator equation on the position bras 
and kets as follows^: 

p,^d{x,y;t) = {x\p,^d{t)\y) =Af{t)exp [~h{t)x'^ - a* {t)y'^ + 2d{t)xy] . (37) 
The master equation (pS)) thus reduces to the following coupled system of differential equations: 

^ = Ahia^-2^{t){ii^ + c) + D{t)-2f{t)~ai (38a) 
^ = 2 [hi - 4 + 5') + ^ H + - 27(i)ai + 2/(t)(aR - 5) (38b) 



dc 
dt 



^aic - 2-f{t){5,ji + c) + D{t) - 2f{t)ai (38c) 



— MM) = 2hi. (38d) 
dt 

Here, the subscript / denotes the imaginary part aj ~ 3[a]. It turns out to be advantageous to directlycompute the 
time evolution of our three non-trivial Gaussian correlators. Analogous methods have been used in [24| to analyse 
decoherence in an upside down simple harmonic oscillator. Let us start by recalling equation (|12p . It is valid for the 
coefficients of the reduced density matrix too. We can thus derive the following set of differential equations: 

^ = -2{ul + n^) - 47(0(|5'> + mt) (39b) 

d{\{x,p)) _ ^^^2 ^ o2^/^2^ , |^2^ / ^ 



= -(^0 + + (r) - f{t) - 27(t) ( ) . (39c 



* Alternatively one could calculate p,-cd using the well known Gaussian propagator of the reduced theory [13 



9 



Initially, we impose that the system is in a pure state: 



'(to)) = ^ (40a) 



2^0 

Wo 

2 



if (to)) = ^ (40b) 



^RP}) = 0. (40c) 

We can then straightforwardly find the phase space area in equation (jISp and the von Neumann entropy for the 
system now follows via equation (fT4|) . We discuss the results after having solved for the statistical propagator in the 
correlator approach to decoherence in the next subsection. 



B. The Correlator Approach to Decoherence 



The statistical propagator in quantum mechanics contains all the relevant information about the three non-trivial 
Gaussian correlators characterising our system as is apparent from equation (jisp . Rather than separately solving for 
these three correlators by means of calculating various expectation values, one can directly solve for the statistical 
propagator. This method is analogous in spirit to the quantum field theoretical computations we performed in 
though there we solved the Kadanoff-Baym equations for the system field to obtain the correlators^. The oscillator 
model of interest is admittedly very simple, so one could use many equivalent methods to solve for the correlators. 
The point of our calculation is that in cases of real physical interest, such as the one examined in it is the only 
tractable method to perform the calculation since for interacting (non-Gaussian) models it is very hard to solve for the 
total density matrix, even perturbatively. It is possible however to perturbatively solve for the statistical propagator. 

As the derivation of the statistical propagator for general initial conditions, i.e.: initial conditions that also permit 
entangled states, is straightforward but rather cumbersome, we present a full derivation in appendix [Xj The principle 
on which this derivation relies is to diagonalise the coupled equations of motion for the two oscillators such that they 
decouple. One can then solve the differential equations trivially, rotate back to the original frame and impose initial 
conditions. It is important to realize that our expression for the statistical propagator in equation (jAlOp is exact. 



1. Initial Conditions I: Pure and Thermal Initial States 



In particular, we are interested in what we henceforth call "pure-thermal" initial conditions. Here, pure refers to the 
initial state of the system, whereas thermal refers to the initial state of the environment field. Since we are dealing with 
only one environmental simple harmonic oscillator, the phrase "thermal initial condition" might be a little awkward. 
We only intend to imply that the environmental statistical propagator probes one point on Planck's distribution. We 
need of course many more oscillators to sensibly talk about a thermal distribution of oscillators. Also, we require that 
the system and environment are separable initially, such that no system-environment correlations are present at tg. 
Using equation (|9b[) . we thus require: 

{iHto)) = 2^ mto)) = f ({i(io),P.(to)}) = 

i'i'ito)) = 2^C0th(^) {plito}) = i^coth(^) {{qito)J,{to)})=0 (41) 

{xito)q{to)) = {pAto)Mto)) = {x{to)pgito)) = {q{to)Mto)) = , 

Having the statistical propagator in equation ()A10|) at our disposal, we can now find the phase space area and the 
resulting entropy in various interesting cases. 



2. Initial Conditions II: Time Translation Invariant States 



It turns out there is a specific class of initial conditions that, as we shall see, do not give rise to oscillatory behaviour 
in e.g. the entropy. Expectation values in this case are a constant despite of the coupling A between the two oscillators. 



^ Indeed such correlator equations also follow on the way to the master equation if one inverts the kernel in the Gaussian exponential of 
the path integral for J in equation lj34|l . The resulting equations are IPI. Usually, they arc obtained by a double Legendrc transform 



Q,[2l,[2i. 
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We present a derivation at the end of appendix [X] Initially, we require: 

(^'(*o)> = 2^ {Plito)) - ^ mto),pAto)}) = 

(9'(*o)) = 2U^ (P^(M) = ^ mto),Mto)})=0 (42) 

{3:{to)q{h)) = (p.(to)p,(io)> = 2^ (•i(io)p,(io)> = (<z(<o)Px(<o)> = , 

Clearly, the system is in a pure state and the environmental oscillator is in a "thermal" state whose temperature is 
fixed by requiring: 

Finally, we see that the initial state is entangled, i.e.: the initial state is not of the form p(0) = Ps{0) ® /5_e(0). 



C. Results 



Let us begin by examining figures [T] and [5J In these figures we show the phase space area as a function of time 
for both the master equation approach to decoherence (gray) and our correlator approach to decoherence (black) . In 
figure [T] the environment is, effectively, at T = as we use f3ujQ = 2000. In figure [2] we use /3a;o — 0.2 such that the 
environmental oscillator is in "thermal equilibrium" at some finite temperature ksT = 5wo- The difference between 
the two phase space areas is pcrturbatively small in both cases (the effective frequency that enters the master equation 
has been calculated by pcrturbativc methods). In other words: both methods agree well. We can thus conclude that 
the entropy one obtains from the perturbative master equation in the Gaussian case in fact stems from neglecting the 
correlation entropy, just as in our correlator approach to decoherence. 

As energy is conserved in our model, the Poincare recurrence theorem applies. This theorem states that our system 
will after a sufficiently long time return to a state arbitrary close to its initial state. The Poincare recurrence time 
is the amount of time this takes. In figure [T] we used almost commensurate eigenfrequencies. Consequently, the 
behaviour of the phase space area as a function of time is very regular. Introducing a non-zero temperature in 
figure [2] automatically makes the eigenfrequencies non-commensurate and thus induces a greater irregularity in the 
time evolution of the phase space area. In other words: the Poincare recurrence time has increased. 

Let us recall equation (fT7|) : 

5'vN ~ S'total = Ss{t) + SE{t) + SsE{t) ■ 

In figures [3l |4] and [5] we show the Gaussian von Neumann entropy for the system Ss{t) (solid black) and for the 
environment Ssit) (dashed black), the Gaussian correlation entropy SsE{t) (dashed gray) and, finally, the entropy 
resulting from the master equation Sl'^{t) (solid gray). The solid black line in figure [3] is the entropy Ss{t) one obtains 
from figure [TJ Likewise, the solid black lines in figures 2] and [5] show the entropy 5*5 (t) as a function of time one 
gets from figure [2] In figure [5] we show the behaviour up to very large times to nicely illustrate the quasi periodicity 
resulting from Poincarc's recurrence theorem. 

The dashed black lines in figures H] and [5] show the environmental entropy Ssit). Just as Ss{t), one can obtain 
Ssit) from the three non-trivial Gaussian correlators (g^), (p^) and l/2{{q,pg}) and by making use of equation (|15p . 
The environmental entropy in figure [3] is precisely equal to the system entropy Ss{t) because both the system and the 
environment are at T = 0. The dashed gray lines in figures [31 H] and [S] show the Gaussian correlation entropy Sssit) 
one obtains from equation (jl9p . Finally, we can easily compare with the reduced von Neumann entropy one obtains 
from solving the master equation (in solid gray). Just as for the phase space area, the results differ only due to a 
small perturbative error in the master equation. 

The oscillatory behaviour in the entropy can be interpreted in two ways: 

• If one averages out the oscillations, the entropy of the system oscillator evolves from 5s (0) = to some small 
non-zero value Ss > 0, where Ss is the (classical) time average Ss = (5s (t)). This average non-zero value of 
the system entropy 5s exists by virtue of a negative correlation entropy. The amount of decoherence that the 
system has experienced thus equals 5s; 

• The Poincare recurrence theorem enforces that the system returns arbitrarily close to its initial state after some 
finite Poincare recurrence time. In quantum mechanics this recurrence time is rather small. This is not the case 
in field theoretical models as we discuss shortly. 
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Figure 1: Phase space area as a function of time. 
Solid black: correlator approach to decoherence. Solid 
gray: master equation approach to decoherence. We use 
uji/ojo = 2, A/oj,^ = 1/2 and /3ujo = 2000. 
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Figure 2: Phase space area as a function of time. 
Solid black: correlator approach to decoherence. Solid 
gray: master equation approach to decoherence. We use 
LJi/ujo ~ 2, A/(jJq — 1/2 and /3a;o = 0.2. 




Figure 3: Entropy as a function of time. Solid black: 
correlator approach to decoherence, which follows from 
equation {T^i. Solid gray: master equation approach to 
decoherence. Dashed gray: correlation entropy, which fol- 
lows from equation it We use uji/uq — 2, A/cjq = 1/2 
and Ploo = 2000. 
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Figure 4-- Entropy as a function of time. Solid black: 
correlator approach to decoherence. Solid gray: master 
equation approach to decoherence. Dashed black: entropy 
of the environment oscillator. Dashed gray: correlation 
entropy. We use lji/loq = 2, X/luq — 1/2 and /3cjo = 0.2. 
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Figure 5: Entropy as a function of time. Solid black: 
correlator approach to decoherence. Solid gray: master 
equation approach to decoherence. Dashed black: entropy 
of the environment oscillator. Dashed gray: correlation 
entropy. We use ui/ujo = 2, \/loq = 1/2 and j3ujQ — 
0.2. The quasi periodic behaviour is due to Poincare's 
recurrence theorem. 



Figure 6: Entropy as a function of time. Solid black: cor- 
relator approach to decoherence. Solid gray: master equa- 
tion approach to decoherence. We use lji/ujo = 201/200, 
X/loq = 1/4 and j3ujo = 0.2 such that we are in the res- 
onant regime. The entropy as obtained from the master 
equation shows unphysical behaviour. 
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Figure 7: Entropy as a function of time. Solid black: en- 
tropy of the system. Dashed black: entropy of the environ- 
ment oscillator. We use uji/ujo = 2, X/ujq = 1/4 and /3 is 
thus fixed via relation The system entropy vanishes 

(pure state) and is time independent. The environmental 
entropy is a constant too. Despite that A 7^ 0, there are 
no oscillations in the expectation values as shown above. 



So far, it is clear that the phase space area or the entropy in our correlator approach to decoherence and in the 
master equation approach to decoherence agree well up to perturbative corrections. Let us now examine the so-called 
resonant regime, where wi ~ in figure IHl Here, we set lji/ujq = 201/200. Our coupling is still small in this case 
so we are well in the perturbative regime where the master equation should yield sensible behaviour. Clearly, it does 
not as the entropy blows up. In the resonant regime the master equation suffers from secular growth which physically 
is not acceptable. Our approach indeed yields a perfectly finite evolution for the entropy. 

In figure [7] wc show the time evolution of the entropy for the time translation invariant initial conditions in equation 
p2)) in our correlator approach to decoherence. Neither the entropy nor the phase space area change in time. We find 
this result rather interesting. Despite the fact that we are dealing with a non-zero coupling A 7^ 0, we can fine tune 
our initial conditions (|42p such that there appears to be no influence of the coupling whatsoever when the entropy is 
measured. In other words, the expectation values in equation (|42p are time independent and consequently fixed by 
their initial value despite of the coupling. This does not depend on how large the coupling A is. 



We are really interested in the general case of A''-oscillators. Let us consider the original Lagrangian in equation 
(|3]). We take N = 50 throughout this section. In this section, we perform again two calculations: 

• Decoherence in the conventional approach: we can numerically solve for the master equation in this case; 

• Decoherence in our correlator approach: we can solve for the three non-trivial Gaussian system correlators by 
exact numerical methods from which we calculate the entropy. 

Again does our simple action in equation (U) allow for a clean comparison of the two approaches to decoherence. 



The form of the Paz-Zurek master equation itself remains unchanged of course and is given in equation (|28p and 
(P^)) . Compared to the previous section where we studied two oscillators only the precise form of the coefficients in 
equation changes. Therefore, we exploit the same strategy as in the previous section and solve for the three 
non-trivial Gaussian correlators as in equation (|39p . 



V. 



N COUPLED OSCILLATORS 



A. The Master Equation Approach to Decoherence 
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B. The Correlator Approach to Decoherence 

The simplest way to study iV = 50 oscillators is not to diagonalise the resulting equation of motion (although this 
is of course possible in principle) , but just study the unitary evolution of all Gaussian correlators in our system of 
interest by numerical methods. Let us recall the total Hamiltonian of our system of interest given in equation 



. „ . . 1 1 ^ 
H^Hs + He + Hint = ^pI + 2^0^^ + 



n— 1 ^ ^ 



Hamilton's equations of motion thus yield: 



X 


= Px 


(44a) 




N 




Px 


= -UjIx - Ang„ 


(44b) 


q-a 




(44c) 


Pq^ 


= -uj^qn - A„i . 


(44d) 



We can thus straightforwardly derive the following set of coupled linear first order differential equations that govern 
the unitary time evolution of all Gaussian correlators: 

^^^'^ - 2(\{i,p^}^ (45a) 
/I \ ^ 

-2loI I ) - 2 ^ XniPxqn) (45b) 



At \2 



At 

'^{\{x,Px}) 



Tl=l 

N 



At 

ml) ,/i 



At \2 
d(pL> 



'^o(a;') + (p')-^A„(xg„) (45c) 

2^^{'Zn,P,„}) (45d) 

-2ujI (^{qn,Pqj'^ - 2Xn{Pq„x) (45e) 

-cjfM) + {pD - >^n{xqn) (45f) 

{Pxqn) + {xpqj (45g) 
AT 

{PxPq„ ) - 1^0 {xqn) -^Xi imn) (45h) 
1=1 

"^^""^'"^ - {PxPqJ - ^lixqn) ~ Xn{x') (45i) 



At 

d(^{g„,Pg„}) 

At 

A{xq„) 
At 

Ajpxqn) 
At 



At 

djPxPqr, 

At 

d(g,ig»i 
At 

At 

d(P<;,.P<;„. 

At 



'oixPqJ - i^liPxqn) - A„ (^{x,px}\ - A„ (|{g„,P,,j\ - ^Xz{qiPq„) (45j) 



1=1 

{QnPqJ + {qmPqJ (45k) 

(P9„P<?,„ ) - ^^m {qnqm) ~ A,„ (xg„) (451) 

(9nj5g„ ) - {q-mPq^ ) - A„ {xpq^ ) - Xm {xpq„ ) ■ (45m) 
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In the last three equations m ^ n is imphed. In other words, for m = n, wc would find equations (j45dp . (|45e|) and 
(j45fj) again. As in the previous section, we require "pure-thermal" initial conditions. They generalise to: 

(i'ito)) = ^ m,)) = ^ {{xihiuto)}) = 

(fnito)) = {plito)) = ^COth(^) ({g„(io),P,„(to)}) =0 (4g) 

(i(io)g„(io)) = {p,{to)pgAto)) = {x{to)pg^{to)) = {qn{to)Px{to)) = 

{qnito)qm{ta)) = {Pq„{io)Pqm{io)) = {qn{to)pq^{tn)) =0, 

where of course 1 < {n, m} < N , m ^ n. For simplicity we will restrict ourselves to the case where A = A„, 1 < n < A^, 
throughout the paper. 



C. Results 



We consider = 50 environmental oscillators. Let us firstly present a case where the master equation produces 
accurate results. In figure [5] wc show the evolution of the entropy obtained from the master equation (solid gray) and 
in our correlator approach to decohcrence (solid black). Wc randomly generate the environmental frequencies a;„/ajo, 
with 1 < ri < A^, in the interval between 2 and 3. Henceforth, we denote this by a;„/ajo e [2,3]. With A/wq ~ 1/8 
this case is clearly deep in the perturbative regime away from resonance frequencies as all environmental frequencies 
are larger than the system frequency. Only a little entropy is generated so virtually no decoherence has taken place 
compared to the thermal value the entropy can reach in principle if decoherence is effective. Indeed, given some 
temperature of the environment, the thermal value of the entropy provides us with a reasonable estimate for the 
upper limit the entropy can reach. For the parameters under consideration, we have a thermal entropy St\i — 0.458 
as I3luq ~ 2, whereas the average generated entropy is only of the order of Ss — 0.08. This case is thus not very 
different from a T = evolution. Our system has not decohered as the environmental frequencies are all larger than 
the system frequency, the coupling is weak and the temperature of the environment is low. In other words, the system 
and environmental oscillators are virtually decoupled from each other. 

Let us now consider figure Here, we used w„/wo G [9/10,11/10] with 1 < n < A^. These frequencies are thus 
randomly spaced around cjo in a small interval. As all frequencies are close to each other, we expect that thermalisation 
and decoherence occur swiftly and effectively, unlike in the previous case we considered. At early times < loqI < 15, 
we see that the entropy calculated from the master equation and our correlator approach to decoherence coincide 
nicely. This is to be expected as the perturbative master equation should yield the correct early time evolution. In 
other words: it is possible to extract the relevant decoherence time scale correctly from both the master equation and 
from our correlator approach to decoherence. 

As entropy or the phase space area quantifies the amount of decoherence, the rate of change of the phase space 
area quantifies the decoherence rate F. We can determine F at each moment in time formally by solving: 

A(t)+F(t)A(t) = 0. (47) 

Due to the oscillatory nature of the phase space area as a function of time, one should average over a conveniently 
chosen time interval to sensibly determine for example the decoherence rate at early times. In this quantum mechanical 
example it is not possible to find a simple (time independent) analytic expression for F, whereas in the field theoretical 
model we investigated F is very well approximated by the tree level decay rate [l| . In [13] one can find a decoherence 
time for this model that depends on the spatial coordinates of the density matrix. We do not find an invariant measure 
of the decoherence rate of that type. 

Now, let us turn our attention to the late time behaviour in figure [5] Around Lu^t ~ 280 the master equation 
suddenly destabilises and the entropy becomes ill defined in the sense that the effective phase space area of the state 
becomes negative. One can easily check that the 51 eigenfrequencies one obtains by rotating to the diagonal frame, 
analogous to equation (|A4[) for A" = 1, are all positive so the resulting evolution should be stable. Clearly, this signals 
a breakdown of the perturbative master equation. 

The latter seems to be a generic feature of the perturbative master equation in the resonant regime. Let us consider 
the figures [TUl [TTl [T^ and [T^ In figure [TU] we use Wn/wo = 1 + n/50, in figure [TT] we use a;„/wo = 1 + n/100, in figure 
[T^ ujn/ujQ G [3/4,3/2] and in figure [T51 ujn/ojo G [19/20,21/20]. Independently on how precisely one chooses the 
distribution of environmental frequencies, the perturbative master equation breaks down. In all these cases arc the 
eigenfrequencies positive, and wc are in all cases in the perturbative regime, so the evolution of the entropy should 
be regular at all times. 

The system's entropy should asymptote to the thermal entropy if thermalisation is complete which we indicate by the 
dashed black line in all the figures in this section. A perfectly thermalised state corresponds to a maximally decohered 
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Figure 8: Entropy as a function of time. Solid black: 
correlator approach to decoherence. Solid gray: master 
equation approach to decoherence. We use a;„/ajo £ [2, 3], 
\/ljq = 1/8 and Pujq = 2. 



s 




100 200 300 400^°^ 



Figure 9: Entropy as a function of time. Solid 
black: correlator approach to decoherence. Solid gray: 
master equation approach to decoherence. We use 
ujn/ojo e [9/10, 11/10], = 1/40 and Pujq = 1- The 

dashed black line indicates complete thermalisation (per- 
fect decoherence). 
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Figure 10: Entropy as a function of time. Solid 
black: correlator approach to decoherence. Solid gray: 
master equation approach to decoherence. We use 
lu„/luo = 1 + n/50, X/loq = 3/40 and /Suo = 2. 
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Figure 11: Entropy as a function of time. Solid 
black: correlator approach to decoherence. Solid gray: 
master equation approach to decoherence. We use 
ijJn/ijJo = 1 + n/100, \/loI = 3/40 and Pu)o = 2. 




Figure 12: Entropy as a function of time. Solid 
black: correlator approach to decoherence. Solid gray: 
master equation approach to decoherence. We use 
lo^/loq e [3/4, 3/2], A/tjg = 1/16 and /3a;o = 2. 




Figure 13: Entropy as a function of time. Solid 
black: correlator approach to decoherence. Solid gray: 
master equation approach to decoherence. We use 
u)u/u)o G [19/20,21/20], A/cjg = 1/10 and fiujo = 1/10. 
At such a high environmental temperature we generate a 
significant amount of entropy. 
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state and consequently an imperfectly thermalised state corresponds to a partially decohered state. Averaging out 
the oscillations of our unitary evolution in e.g. figure [TU] yields an average value of the entropy that is below this 
thermal value. This can be expected given the fact that our coupling is very small A/cjq ~ 3/40. Also, one can verify 
that the energy is conserved by the evolution as it should. 

However, the master equation in figure[TO]yields a stationary entropy that is larger than the thermal value, indicating 
that the temperature of the system would be higher than the environmental temperature. Clearly, this does not make 
sense. The two most important quantitative measures of decohcrence one hopes to extract out of an experiment or 
a calculation are the decoherence rate and the total generated entropy as the latter tells us how much decoherence 
has taken place and how classical the state has become, and the former how fast this state is reached. We conclude 
that for this choice of parameters the total generated entropy at late times does not follow correctly from the master 
equation. We checked that this failure of the master equation is generic. Even in the case of figure IHl we see that 
although the asymptote is roughly equal to the thermal value in the range 25 < wot < 280, it is too high based on 
our exact numerical analysis {Ss < Sth). 

The oscillations in the entropy as a function of time that are present in our exact evolution are not present in 
the evolution of the entropy that follows from the master equation. The master equation thus tends to overdamp 
oscillations in expectation values. 

Finally, compared to the iV = 1 case we previously considered, we observe that Poincare's recurrence time has 
dramatically increased. For example in figure [T^] one has to wait for a very long time before a random fiuctuation 
decreases the entropy significantly again to values close to its initial value. Thus, by including more oscillators 
decoherence becomes rapidly more irreversible, as one would expect. 

VI. CONCLUSION 

Decoherence is often studied by considering a perturbative master equation. This approach suffers from both 
theoretical and practical shortcomings. It is unsatisfactory that pred evolves non-unitarily while the underlying theory 
(quantum mechanics or quantum field theory) is unitary. We are not against non-unitary equations or approximations 
in principle, however, one should make sure that the essential physical features of the system one is describing are 
kept. The master equations does not break unitarity correctly, as we have shown in this paper. On the practical side, 
the master equation is so complex that field theoretical questions have barely been addressed: there does not exist 
a treatment to take perturbative interactions properly into account, nor has any reduced density matrix ever been 
renormaliscd. Moreover, as we show in this paper, the perturbative master equations leads even in simple situations 
in quantum mechanics to physically unacceptable behaviour. Although for early times the entropy increases, for 
late times the approach fails and suffers from secular growth as the entropy continues to grow without bound. This 
behaviour is generic when at least some of the environmental oscillators are in the resonant regime. Of course one 
could argue that we should have solved the non-perturbative master equation in the resonant regime. Although true 
in principle, the point of this paper is that it is much more convenient on theoretical and practical grounds to solve 
for correlators than for density matrices. 

We advocate a novel approach to decoherence based on neglecting information in observationally unaccessible 
correlators. If proper interactions (non-Gaussianities) are taken into account, neglecting the information stored in 
higher order no n- Gaussian correlators gives rise to an increase in the Gaussian von Neumann entropy. In this paper 
we consider a quadratic model of N coupled simple harmonic oscillators where no non-Gaussianities are generated. 
Due to the non-zero coupling however, an increase in the Gaussian von Neumann entropy can be observed at the 
expense of a negative correlation entropy as the total von Neumann entropy is constant in unitary theories. 

In this paper we study a quadratic model where the trace can be performed, where the perturbative master equation 
can be solved and where there is no need to renormalise as we study decoherence in a quantum mechanical setting. We 
can thus circumvent the practical drawbacks that usually prevent us from solving the perturbative master equation in 
a renormaliscd interacting quantum field theoretical model. In this simple quadratic quantum mechanical model, we 
can thus actually compare the master equation approach to decoherence and our correlator approach to decoherence. 
From the numerical analysis we conclude: 

• Away from the resonant regime, where all a;„ significantly differ from loq^ the evolution of the entropy in the 
master equation approach and in our correlator approach to decoherence agree nicely (up to the expected 
perturbative corrections). The system decoheres only weakly in this regime however. The two approaches also 
agree nicely at very early times such that both approaches can be used to accurately calculate decoherence rates; 

• In the resonant regime, where a certain number of aj„ ~ loq^ the perturbative master equation breaks down 
despite of the fact that we are still in the perturbative regime and despite of the fact that all eigenfrequencics 
are positive. The value to which the entropy asymptotes in the master equation approach before the breakdown 
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is reached is generically too high. This leads to an inaccurate prediction of the amount of decoherence that has 
taken place. In some cases, the breakdown occurs already before this asymptotic state has been reached; 

• The evolution of the entropy in both the resonant and in the non-resonant regime in our correlator approach to 
decoherence behaves perfectly finite. The late time asymptote that is eventually reached either indicates perfect 
thermalisation (maximal decoherence) or, for smaller values of the coupling and a finite number of environmental 
oscillators, imperfect thermalisation (imperfect decoherence). 

Our correlator approach to decoherence thus provides us with a new insight in the conventional approach to de- 
coherence. If the entropy from the master equation agrees with our exact Gaussian von Neumann entropy up to 
perturbative corrections, then one neglects the information stored in the system-environment correlators in the con- 
ventional approach to decoherence too. 

There are generally speaking two quantitative results that can be obtained after performing an experiment or a 
calculation concerning decoherence: 

• Decoherence rate; 

• Total amount of decoherence. 

We argue that one should use the entropy (or, equivalently, the phase space area) to quantify these two aspects of 
decoherence. The decoherence rate F is the rate at which the phase space area changes and can formally be obtained 
from solving equation (|47[) . The total amount of entropy generated measures the total amount of decoherence that 
has occurred. At late times we have seen that, upon averaging out the oscillations, a thermal state is maximally 
decohered. This agrees perfectly with the point of view where a reduced density matrix whose off diagonal elements 
(in the observer's basis) have disappeared corresponds to a classical state. A (free) thermal density matrix follows 
from equations (|10bp and (fT3)) as: 



pthix,y;t) =7Vth(i)cxp 



(48) 



The thermal correlators can be read off for example from equation (|4ip . Clearly, at high temperatures we have 
that Ath ^ 1 such that the off diagonal terms in the density matrix have almost disappeared. The classical limit of 
quantum mechanics in not classical mechanics but rather classical stochastic mechanics. Likewise the classical limit of 
quantum field theory (in a thermal environment) is not classical field theory but classical stochastic field theory. This 
implies that a particular measurement does not yield a certain predetermined outcome but rather a certain outcome 
that is randomly drawn from a probability distribution function of uncorrelated or not entangled possibilities^. 

We can learn several things about decoherence in field theories from our quantum mechanical analysis. Firstly, 
the reason why the perturbative master equation fails in the resonant regime is that the quadratic coupling between 
system and environment is treated on equal footing as a proper non-Gaussian interaction. In the perturbative master 
equation, one attempts to solve the following "interaction" in the "Kadanoff-Baym" equations: 

+ + . . . . (49) 



Of course, this is not a proper interaction and the 1 particle irreducible (IPI) self-mass that enters the Kadanoff- 
Baym equations'' due to such a quadratic coupling vanishes. In this simple quantum mechanical example, one 
should therefore just solve for the correlators following from the von Neumann equation exactly as we showed in 
this paper. We thus also accounted for the backreaction from the system on the environment. The set of Kadanoff- 
Baym equations is the sophisticated machinery to properly solve for the propagators in an interacting quantum field 
theory. Only the diagrams that follow from a 2PI effective action contribute. As is well known, the Kadanoff-Baym 
equations resum Feynman diagrams and prevent secular terms from developing. This leads to a stable, thermalised 
late time behaviour (see 0, [H, H^l)- Secondly, for N = 1 and for almost commensurate eigenfrequencies, we have 
seen that Poincare's recurrence time is rather small: the system returns rather quickly to its initial state. For 
non-commensurate eigenfrequencies Poincare's recurrence time increases. For = 50 we have seen that Poincare's 



The Wigner transform of a density matrix is essentially a Fourier transform with respect to the relative coordinate x — y in e.g. equation 
II48I I such that the cross section of the Wigner transform of the thermal density matrix is a circle. Intuitively, Wigner space provides us 
(almost) with a classical stochastic probability distribution function on phase space. 

The Kadanoff-Baym equations are the equations of motion for the various propagators in the in-in formalism in quantum field theory 
that stem from a 2PI effective action such that the contributing diagrams are IPI. 
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recurrence time dramatically increases such that we have to wait for a much longer time before the system returns 
to a state arbitrarily close to its initial state. There are two complementary ways to interpret this observation. One 
could take the point of view that our system has experienced an average entropy increase, e.g. < 6*5 ^ S'th, where 
Ss is the classical time average of Ss{t). In other words, the system has experienced a total amount of decoherence 
equal to Ss- Alternatively, one can say that since the Poincare recurrence time is still finite, although large, no 
irreversible process of decoherence has taken place as the system returns arbitrarily close to its initial state again. In 
an interacting field theory several modes couple due to the loop integrals (hence iV — >■ 00) and clearly our Poincare 
recurrence time becomes infinite. Hence, the entropy increase has become irreversible (for all practical purposes) and 
our system has decohered. 

We hope to generalise the analysis presented here to a quantum mechanical interaction that comes closer to field 
theory but is still free from divergences. In particular, we will consider entropy generation in a quantum mechanical 
Xxq^ model, which, when analysed by making use of the master equation, is very similar to the present Xxq model 
after q has been integrated out [ll|. The purpose of introducing our novel approach to decoherence is not to study 
quadratic quantum mechanical models, but proper interacting field theories, like in [l|. 



Appendix A: The Statistical Propagator for A'^ = 1 



Let us now return to the lagrangian (|36|) from which it follows that the equation of motion for the two oscillators 
can be written in matrix form as: 

^q + nq^O, (Al) 



n=[i^ ) (A2) 



where: 

Al ujf 

and where we employed the notation q"^ — {x,q). We can now easily diagonalise the equation of motion: 
where we diagonalised the system by making use of the rotation matrix R 



+ ^dqd = , (A3) 



* wry V i(c^2^C^2)^ l^(^2_^2)2+4^2 / 



2 



\ /i(w§+^?)-lVR^^iF+4Af 



with: 



Moreover, we defined: 



Finally, one can derive that: 



R^[ -^";(^) ) . (A5) 

' sm(6') cos(t^) ' ^ ' 



r7 i 1 - f cos(6l)x(i) - sin{0)q{t) » , , 

'^'^^^'^ I q{t) J ^ \sm{e)x{t)+cos{e)q{t) ' ' 



cos{2e) = , ^ " = (A7a) 

sin(2^) = I =. (A7b) 

In order for the system to be stable we should have ujq > 0, which implies Ai < loquji. When Ai > ujquji the upper 
eigenmode corresponds to an inverted harmonic oscillator, ujq < 0, and we shall not consider this case here. We can 
immediately solve the equation of motion in the diagonal frame: 

rr (A ^ f ^(*) \ - f ^0 cos(a)oi) + Bosm{u}ot) \ , , 

^'^^^ ~ I qit) J cos(wit) + Bi sin(LSit) ) ' ^^^> 
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We do not yet impose initial conditions for the coefficients in this solution, but rather first solve for the statistical 
propagator for the system: 



F4t;t') = ^Tr[p{to){x{t'),m}]^l{m'),m}) 



(A9) 



cos'{e){{S:{t'),S:{t)})+sm\e){{q{t'),m})+cos{e)sm{0)a^^^^ 



Inserting the solution (|A8|) yields the general form for the statistical propagator: 



(AlO) 



cos' 



'(6')|(Ao) cos(woO cos(woi') + (Bq) sin(cLiot) sin(a;oi') 



{{Ao,Bo}) 



{cos{ujot) sm{iOQt') + sin(u;oO cos(wot')) | 



+ sin^(6')|(ii) cos(£:jii) cos(wit') + (Bf) sin(wit) sin{u)it') + ^i^i^-^i^) (cos((:jit) sin(a;it') + sin(wii) cos((I)it'))| 
^ sin^20) 1^1^^^ Ai}) (cos(wot) cos(LJit') + cos((Dit) cos(woi')) + ({-^o, ^i}) (sin(ajot) sm{uiit') + sin(ajit) sin((I)ot')) 
+ ({Ao, Bi}) (cos(ajoi) sin{ujit') + sin(ajit) cos(wot')) + ({-Bo, ^i}) (sin((Dot) cos{LUit') + cos((Dii) sin((Dot')) |. 



Similarly, we can derive the statistical propagator for the environment: 



1 



1 



(All) 



sin2(0)({i(t'), m}) + coa\e){{qit'), m}) - cos(0) sin(0)({l(t'), 1(0} + {!(*'), ^(01) 



such that we find: 



(A12) 



3in^(6')|(Ao) cos(a;oi) cos(woi') + (Bq) sin(a;oi) sin(iuoi') + 

(6')|(A^) cos(ajii)cos(ajit') + (i?J) sin(wit) sin(wi<') + 



{{Ao,Bo}) 



' cos 
sin(26') 



{{Ai,B,}) 



{cos{u!Qt) sin(a;oi') + sin(a;ot) cos(a;ot')) | 
(cos(wit) si-a{ujit') + sin(a;it) cos(ajit'))| 



|({Ao,^i}) (cos(woO cos(ajii') + cos(wi<) cos(woi')) + {{^o,Bi}) (sin(woO sin(wit') + sin(a;it) sin(cjot')) 
+ ({Ao, Bi}) (cos(ajoi) sin(cJii') + sin{uiit) cos(a)oi')) + ({-^o, ^i}) (sin(cJoO cos(wii') + cos(a;i<) sin(a)oi')) |- 



Finally, we need the statistical propagator for the system-environment correlations: 



1 



1 



(A13) 



sin(26l) 



{{{m,m}) {{Ht'),Ht)})) +co,^9){m'),m}) ~^in'io){m'),m) 



which in turn yields 
sm{26) 



(A14) 



I (Al) cos(a;it) cos(wit') + {Bf) sin(a;i<) sin(wit') + 



{{AuB,}) 



(cos(wit) sin(a;it') + sin(a;ii) cos(wit')) 



'{Al) cos(a;ot) cos(woi') — {Bq) sin(a;o<) sin(wot') — 



{{Ao,Bo}) 



(cos(woO sin(a;oi') + sin(wot) cos(cjot'))| 



cos'^(6) r . - . . 

+ — ^|({Ao, Al}) cos(wit) cos(a)oi') + ({^o, Bi}) sin((I;it) sin(woi') 

+ {{Ao, Bi}) sm{u}it) cos(cjoi') + ({So, A}) cos(wi<) sin(wot')} 

({io, Al}) cos(woi) cos(£:jit') + ({Bo, Bi}) sin(cjot) sin(wit') 

+ ({Ao, Bi}) cos(woi) sin(a'ii') + ({Bq, ^i}) sin(woO cos 
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and Fqx{t;t') = Fxq{t']t). Now, we can impose initial conditions at on all initial correlations. We can derive: 



(i'(io)> 
('z'(io)) 
{x{h)q{h)) 
iPlito)) 

{Pxito)Pq{to)) 

({x(to),p,(to)}> 
({g(to),Pg(to)}) 

{x{to)Pq{tQ)) 

{q{to)px{tQ)) 



cos\e){Al)+sin'{9){Al) 



sin' (0) (A'a) +cos'i0){ Al 

1 
2 

2 / /) V— 2 / D 2 \ I . . n n 2 / /) \ — 2 / r> 2 



sin(20)({io,ii}) 
sin(20)({io,ii}) 



sin(20) ((yl?) - {Al 



1 
2 

cos^ 



cos(20)({Ao,Ai}) 

1 

2 
1 

2 



sin(20)woWi({So,Bi}) 
sin(20)cSo<:^i({Bo,Bi}) 
cos{20)ojoU'i{{Bo,Bi}) 



(Al5a) 
(A15b) 
(Al5c) 
(A15d) 
(A15e) 
(A15f) 



sin(20) (^ojf{Bf)-ojUB'o) 
i0)Qo{{Ao,Bo})+siu'i0)ui{{Ai,Bi}) + ^sm{20) (cSi({io, Bj}) + ^o({^o, ii})) (A15^ 



sm'{0)Ldo{{Ao,Bo})+cos'{e)Qi{{AuBi}) - - sin(20) c^i ({Ao, Bi}) + cDo({Bo, Ai}) (A15h) 



sin(26l) 
2 

sin(26') 



u;i{{Ai,Bi}) - LUo{{Ao, Bo})j + CDs'" {0)lui{{Ao, B^}) - sin^{9)Cdo{{Bo, ii}jjA15i) 
Cji{{Ai,Bi}) - wo({io, Bo})) + cos2(0)cSo({So, ii}) - sin2((?)cDi({io, BilfAlSj) 



These equations can be inverted to give: 
1 



{AD 

({^,^}) 

(B'o) 
0!) 
{{Bo,B,}) 

{{Aq,Bo}) 



2 ^.x'ito)) + (q^to)) + cos(20) {{x'ito)) - if (to))) -2sM20){xito)qito)) 

i [{x'{to)) + (q^ito)) - cos(20) ((^2(to)) - if (to))) + 2sin(20)(£(to)g(to)> 
- sin(20) + 2 cos(20)(a;(io)g(io)) 



1 

1 

2aSf 
1 



(P'(io)) + (p2(to)> +cos(2^^) {{Plito)) - {Plito))) -2sm{20){pMP,ito)) 
\pl{t^)) + {f,{to)) - cos(20) {{plito)) - ifqito))) +2sm{20){p,{t,)pq{t,)) 



UJQUJI 

1 

2uio 



({Ai,Bi}) = — 



1 

2^1 



{{Ao,B,}) = ^ 



({^o,ii}> = 



1 

2^0 



((p2(to)) - ifqito))) sin(20) +2cos(20)(p,(to)p,(to)> 
({i(to),P.(to)}> + ({g(io),?5,(to)}) + cos(20) {i{x{to),Mto)}) - i{qito),Pg{to)})) 
-2sin(20) iiqito)pAto)) + ixito)pqito))) 

i{xito),pM}) + i{q{to),Pgito)}) - COS(20) ii{x{to),Mto)}) - i{qito),Pg{to)})) 

f2sin(20) iiqito)Pxito)) + ixito)Pqito))) 

sin(20) ii{xito),Mto)}) - i{qito),Pqito)})) + 2 ((x(io)p,(<o)) - (g(to)p..(io))) 
+2cos(20) iixito)pqito)) + iq{to)Pxito))) 

sin(20) ii{xito),Mto)}) - i{qito),Pgito)})) -2iixito)pqito)) - iqito)Mto))) 
h2cos(20) {ix{to)pq{to)) + iq{to)pAto))) 



(Al6a) 

(A16b) 
(Al6c) 
(A16d) 

(A16e) 

(A16f) 

(A16g) 
(A16h) 
(Al6i) 
(Al6j) 



Note that for example x and q commute. It might seem that we brought out the big guns to solve such a simple 
problem. This is necessary, however, to generalise the standard setup to include non-separable initial states. Both Paz 
and Zurek and Caldeira and Leggett [l^ for example assume that the density matrix is separable separability. 
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p(0) = ps(0)®/5£(0). 



(A17) 



As is apparent from equations (jAlsp above, we are now in the position to easily relax this assumption and generalise 
to non-scparablc initial states. Grabert et al. [ist and Romero and Paz [28j do not assume separable initial conditions 
in their discussion of the master equation and consider more general initial states too. 

For the "pure-thermal" initial conditions given in equation ()41|) . we obtain by means of equation (jA16|: 



{AD 

{AD 
{{A,.M}) 
{Bl) 

{Bl) 
({So, Si}) 



1 1 

1~ - — coth 

zwq 2a; 1 



2 



cos(26') 



1 



1 



2ujQ 2a; 1 



■ coth 



2 



1 



1 



sin(26i) 



2a;o 2a; 1 



/3uj 



1 








2 




1 


ujo 






2 




1 


- sin 







\ 2ujo 2uji 
coth 

coth 



2 

coth 



- — + - — coth -— ^ - cos(26') — coth -— ^ 



1 



1 



2a;o 2a; 1 



2 



2 



cos(26i) 
(^)-cos(2.)( 



uji 



— coth 



2 2 



coth 



2 
2 



-Z coth 



(A18a) 

(A18b) 
(Al8c) 
(AlSd) 
(Al8e) 
(AlSf) 



with all other correlators vanishing. 

Let us now discuss the time translation invariant states as initial conditions as discussed in subsection IIV B 21 In 
order to investigate the dependence on the initial conditions, it turns out to be advantageous to rewrite equation 
(jAlOp in terms of the average time and the time difference: 



1 



{t + t') 



At = t-t' . 

Making use of several trigonometric identities, this yields: 



(Al9a) 
(A19b) 



F^ir-At) = 
cos2(6') 



2 

,2/ 



sin^(6') 
2 

sin(26l) 



cos(LSoAt) [(ig) + {Bl)\ +cos(2wot) [(ig) - (B^)J + sin(2woT)({io, Bo}) 
cos(a;iAt) [(i?) + + cos(2a;ir) [(i?) _ + sin(2a}ir)({ii, Bi}) 



(A20) 



cos{Aujt) cos(wAO [{{Ao, Al}) + ({Bq, Si}) j + cos(2wr) cos(AwAi/2) [{{Aq,Ai}) - {{Bq, BJ) 
sin(Awr) cos(wAt) (({Bo,^}) - ({ic^i})) + sin(2wr) cos(AwAV2) (({-Bo,ii}) + ({io,Si}) 



where we defined: 



1 



LU = -(a;o +a;i) 



(A21a) 

Aw dio - wi . (A21b) 
We now require that the statistical propagator of our system does not depend on the average time r, such that: 



{AD 
{A\) 



{B', 



(A22a) 
(A22b) 
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and all other correlators should vanish. Equation (|A15[) thus tells us: 



(£'(to)) 




(A23a) 


ifito)) 




(A23b) 


{x{to)q{to)} 


= -sin{29)[{Al)^{Al)) 


(A23c) 


(Plito)) 


= coa^i0)Ldl{Al) +sm^i0)Ldf{Al) 

\/U\ U/ V / i\ 1/ 


(A23d) 


\Fq \^() J/ 


^ sm^{e)Ql{Al) +cos^i0)Ldl{Al) 


(A23e) 


{Px{to)Pq{to)) 


= -sin{20){u>l{Al)-col{Al)) 


(A23f) 


{{x{to),p^{to)}) 


= 


(A23g) 




^ 


(A23h) 


{x{to)Pq{to)) 


= 


(A23i) 


{q{tn)Px{to)) 


= 0. 


(A23j) 



We have thus found a 2-parameter family of initial conditions such that the statistical propagator for the system does 
not depend on the average time. Consequently, the entropy does not depend on the average time. Let us require that 
our system is in a pure state initially, such that: 



{x'ito)) = (A24a) 

ZLOq 

(Plito)) = ^, (A24b) 



2 

which is equivalent to: 



2wo 

The other correlators in (|A23|) can now trivially be determined: 



{Al) = {Al) = ^. (A25) 



if (to)) = 7^ (A26a) 



2w, 

2 

U! 



'0 



imo)) = TT- (A26b) 

{x{to)q{to)) = (A26c) 

{pMPqih)) - (A26d) 

Clearly, the environment is in a thermal state at a temperature dictated by requiring equation (j43|) : 



coth 



V 2 y t^o 



Appendix B: Reducing the Density Matrix 



In this section we compute the reduced density matrix for our two coupled simple harmonic oscillators. Starting 
point for this derivation is the full density matrix that contains both the system and environmental degrees of freedom 
in the Hamiltonian 



{x,q\p{t)\y,r) = p{x,q;y,r]t) = TVexp 



X q 



+2(.T q 



a di/2 
di/2 ai 

c ei/2 
el/2 ci 



y r 



a* d\/2 
dl/2 a\ 



(Bl) 
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where c and Ci are real such that there are precisely ten degrees of freedom in this density matrix. Of course the 
density matrix is hermitian: p(x, q; y, r; t) = p*(jj, r; x, q; t). The normalisation constant J\f is determined by requiring 
Tr[p(i)] = 1, yielding: 



where, for normalisability we required: 



af-ci > (B3a) 
A{a^ ~c){a^-ci) > {df - eff . (B3b) 



The unitary dynamics of the ten degrees of freedom in this density matrix is governed by the von Neumann equation 
([2]). Of course we can trace the density matrix which yields: 

/oo 
dqp{x,q;y,q;t) :^Ncxp -ax^ -hy"^ + 2cxy , (B4) 
-OO 



where: 



- 8{af - c,)^''' - 


eif 








Mi-ei|2 




8(a? -ci) 





(B5a) 
(B5b) 
(B5c) 



V2(af-ci) Y TT 27r(af-ci) ^ ^ 

The total von Neumann entropy can now straightforwardly be obtained using the replica trick 0] : 



where: 

X 2 



(B7) 



Although this expression is the final answer for the entropy generated by tracing out the environmental degrees of 
freedom, it is not in a convenient form to study the dynamics. All we have done is relate the von Neumann entropy 
to the coefficients in the full density matrix. The strategy is as follows: the ten degrees of freedom in the density 
matrix can straightforwardly obtained by numerically solving the full density matrix. If we insert the Ansatz for the 
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density matrix (|B1|) in the von Neumann equation, we find the following set of differential equations: 



= iaian+idfdl-efel) (B8a) 

^ = ^ + 2 (a? - 4 + - ^ {{dfr - iefr {d\r + {e\f) (B8b) 

^ = 4aic - (df el - ef 4) (B8c) 

^ 4alaf + (df 4 - ef el) (B8d) 

^ = ^ + 2 {{a\f {afr +cl)-l {{dfr {efr {d\r + {em (B8e) 

^ ^ Aa\c, + {dfe\+efd\) (B8f) 

^ = 2 [(aR + af) d\ + (ai + a\) df + {c - c^) e\] (B8g) 

^ = A + 2[-(aR, + af)d« + (ai + al)4 + (c + ci)ef] (B8h) 

^ = 2 [(aR, - af) e\ + (ai + a\) ef + (c + ci) 4] (B8i) 

^ = 2 [- (or - af ) ef + (ai + a\) e\ + {c-c,) df] . (B8j) 



For pure and thermal initial states as previously diseussed (the system is in a pure state, the environment in a thermal 
state), one can straightforwardly show that the reduced Gaussian von Neumann entropy that results from equation 
()B6|) coincides precisely with the Gaussian von Neumann entropy Ss{t) defined in equation ()18p . We can thus confirm 
the identity (PT|) in an explicit model: 

Ss{t) = Slt{t). 
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